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ABSTRACT 

Wc present the first description of mode trapping for sdO models. Mode trapping 
of gravity modes caused by the He/H chemical transition is found for a particular 
model, providing a selection effect for high radial order trapped modes. Low- and 
intermediate-radial order p-modes (mixed modes with a majority of nodes in the P- 
modc region) are found to be trapped by the C-O/Hc transition, but with no significant 
effects on the driving. This region seems to have also a subtle effect on the trapping 
of low radial order g-modes (mixed modes with a majority of nodes in the G-modc 
region), but again with no effect on the driving. We found that for mode trapping to 
have an influence on the driving of sdO modes (1) the mode should be trapped in a 
way that the amplitude of the cigenfunctions is lower in a damping region and (2) in 
this damping region significant energy interchange has to be produced. 

Key words: stars: oscillations - stars: variables: other 



1 INTRODUCTION 

Subdwarf O (sdO) stars cover a domain of approximately 
60 000 K in T e g and 2.5 dex in \ogg in the hot and 
high gravity domain of the HR diagram. Due to the wide 
spread of physical parameters their evolutionary stage is 
not completely understood, but some progress has recently 
been made with larger than ever (c omprising up to a 
hundred sdOs) spec troscopic studies (|Stroer et al.l 120071 : 
iHirsch fc Hebeill2008l ). D etermination s of T e g , log g and he- 
lium abundance allowed IStroer et al. | (|2007l ) to make some 
statistical correlations between position in the HR and evo- 
lutionary stage, and as a result they propose that sdOs 
come from a mixed population of post-EHB (i.e. sdB de- 
scendants), post-AGB and post-RGB objects. The canonical 
picture depicts them as objects of about half a solar mass 
with an inert carbon-oxygen core, and helium and hydrogen 
burning layers. 

A milestone for sdOs was passed in 2006 with the 
discovery of the first, and unique to date, pulsating sdO 
SP SS J160043. 6+07480 2. 9 (known in short as J1600+0748) 
bv lWoudt et all 1200^. This object is a very fast pulsator, 
with a m ain period of ^120 s, and at le ast seven other fre- 
quencies (|Rodrlguez-L6pez et alj|2009al ) between the main 
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frequency and its first harmonic . The excited frequencies 
were identified bv lFontaine et all (|2008l ) as p-modes and its 
pulsations explained as a classical k mechanism with the 
aid of appropriate models including radiative levitation of 
iron. Our own attempts of exciting p- mode pulsations in uni- 
form m etallicity sdO models failed (|Rodriguez-L6pez et al.l 
l2009bl . Paper I), although some models, as this one, were 
found to present a certain tendency to instability in the p- 
mode region. During the course of our investigation of pul- 
sational properties of sdO models, for which g-mode driving 
was found, a particular model presenting g-mode trapping 
was discovered. 

In main sequence stars the asteroseismic potential of 
g-mode trapping to study the mixing processes within the 
convective core and their effects on the surro unding chem- 
ical g radient have been well documented by iMiglio et al.l 
(|2008l ). In the case of compact pulsators, q-mode trapping 
has b een described extensively for sdBs (|Charpinet et al] 
|2000| ). DA and DB white dwarfs, and GW Vir pre- 



white dwarfs (see e.g. IWinget. Van Horn fc Hansen! Il98ll : 
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Brassard et al.lll99ll; iBrassard et al, 1 1 19921 . or more recently 
Corsico fc Altha us 2006; iGautschv fc AlthausH2002T ) as due 



to chemical stratification in the envelope of these stars. Be- 
cause gravity modes are trapped in the outer envelope, they 
have lower amplitude in the damping core, which intensi- 
fies their instability. This may help, in some cases, alleviate 
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Table 1. Main physical parameters and mass fractions of the sdO model used. X(other) refers to the mass fractions of all the other 
elements. Z is the initial and current metallicity of the model. 



T eS log g 


M 


rm X(H) 


X(He 3 ) 


X(He 4 ) 


X(C) 


X(N) 


X(O) 


X(other) Z 


(K) 


(M ) 
















45 000 5.26 


0.469 


0.650 0.59 


1.2E-03 


0.35 


7.4E-03 


4.2E-03 


2.4E-02 


1.5E-02 0.05 



the problem of observing less modes than predicted by the 
models. It also provides a useful tool t o determine the mass 
and width of the surface la yer (see e.g. Brassard et al.|[l992l ; 
iKawaler fc Bradley|[l993 ). lcharpinet et al.l (|200oF reported 
gravity mode trapping due to the transition between the 
helium core and the hydrogen-rich envelope in sdB models. 
They also identified subtle, but non-negligible departures 
from uniform frequency spacing for pressure modes, caused 
by the same chemical transition region. 

In sdB models, the maximum gradient in Brunt- Vaisala 
(BV) frequency is produced at the transition of the helium 
radiative core to the hydrogen envelope. In sdO models the 
BV profile is more complex, since as evolution proceeds, a 
C-O core builds up, resulting in two chemical transitions: 
one from the C-O core to the He burning shell, and another 
one from the He radiative shell to a H burning shell. This 
renders the mode trapping effects more complex in sdO than 
in sdB models. 

We report in this study the discovery of g^mode trap- 
ping in a sdO model caused by the He/H transition. This 
provides a selection mechanism in the way that trapped 
modes would be more easily excited. Low- and intermedi- 
ate radial order p-mode£] are found to be trapped by the 
deeper transition from the C-O/He, but with no significant 
effects on the driving. 

The importance of g-mode trapping is evident as a way 
to probe the deep stellar interior, where gravity modes prop- 
agate, unattainable in any other way; but also to derive 
information about the location and width of the chemical 
transitions and the mass of the envelope. 

The paper is organized as follows: in Section 2 we de- 
scribe some general characteristics of the sdO model used. 
In Section 3 we describe the (/-mode trapping and the effects 
of cancelling out the He/H and C-O/He chemical transitions 
in BV frequency on the trapping and driving of the modes. 
Section 4 does a similar treatment for p-modes. Finally, Sec- 
tion 5 presents a summary and Section 6 the discussion and 
conclusions. 



2 THE MODEL 

This theoretical exercise arose as an exploration of a particu- 
lar sdO model, whose properties were thoroughly described 
in Paper JJ and parameters given in Table [T] The model, 



1 We remind the reader that in sdO models, modes with low to 
intermediate frequencies behave as mixed modes, i.e. have nodes 
both in the P- and G-mode region. The classification 'pressure 
modes' refers then to mixed modes with most of their nodes in 
the P-mode region (see fig. 7 of Paper I). 

2 The model corresponds to model 8, also named eta650t45, ac- 
cording to its mass loss rate parameter and effective temperature. 
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Figure 1. Left: Normalised nuclear energy generation rate as 
a function of the fractional mass depth parameter (log q = 
1 — M r /Mj<). The deeper peak at logq ~ —0.8 corresponds to 
the helium burning shell; the shallower peak at logq ~ —2.2 
to the hydrogen burning shell. Right: Brunt- Vaisala frequency 
(solid line) showing the two peaks corresponding to the maxi- 
mum chemical gradients: from the C-O core to the He burning 
shell (log q ~ —0.8), and from the He radiative shell to hydrogen 
shell burning (logq ~ —2.2). Lamb frequency (dashed line) for 
1 = 2 modes is also shown. 



built with JMSTAR code jLawlor fc MacDonaldlboOot ) . comes 
from a 1 MQ star on the pre-main sequence with enhanced 
mass loss rate on the red giant branch, that drives the star 
to evolve to higher temperatures at constant luminosity. A 
delayed helium-flash puts the star back in the horizontal 
branch, and further evolution brings the star to the sdO do- 
main. The sdO models have developed a carbon-oxygen core 
while helium and hydrogen shell burning is still produced 
(Fig. [1] left). The change in chemical composition from the 
C-O core to the He burning shell, and from the He radiative 
layer to the H burning shell is seen in the two steep peaks 
in BV frequency and subtler transitions in Lamb frequency 
(Fig. m right). 

The model, which was subject to a non adiabatic ana l- 
ysis with graco (|Mova et alj|2004 iMova fc Garridoll2008l ). 
was found stable in the frequency range ~0.5 to 25 mHz. 
However, an oscillatory behaviour of the growth rate at low 
and intermediate frequencies (see fig. 10 of Paper I) caught 
our attention. That drove us to explore the behaviour of the 
kinetic energy of the modes, which revealed that the model 
experienced mode trapping effects. 
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Figure 2. Logarithm of the total kinetic energy vs. radial order 
for g-modes (this and subsequent plots have been made for i = 
2 modes) of the original (upper circles) and perturbed models 
without the He/H (upper solid line) and C-O/He (upper crosses) 
transition regions. Lower symbols plot the corresponding growth 
rates whose values can be read on the right axis. 
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Figure 3. Period spacing for g-modes with the same degree I 
and consecutive radial orders for the original (circles and dotted 
line) and perturbed He/H (solid line) and C-O/He (crosses and 
dashed line) models. 



3 GRAVITY MODES 

The logarithm of the total kinetic energy of each p-mode 
and its growth rate r\ are given in Fig. (upper and lower 
circles, respectively, the meaning of the other symbols will 
be explained below). The plot reveals that in the region of 
very high-radial order g- modes (radial order k > 40), there 
are small variations in the kinetic energy, of no more than 
one order of magnitude, between different modes. We also 
see that modes with local minimum values of the kinetic 
energy have local maximum growth rates and vice versa. 

The non-uniform distribution of the kinetic energy re- 
sembles the mode trapping phenomenon caused by the po- 
tential barrier due to the composition transition regions. The 
chemical transition regions act as resonant cavities for the 



modes, changing the amplitudes of the eigenfunctions in dif- 
ferent regions of the star and hence modifying their kinetic 
energy, which is given by: 



F 1 n* 



(I6-I 2 + 1(1 + 1) \i h \ 2 )^pr 2 dr 



(1) 



where a is the eigenfrequency, £ r and £j, the radial and hor- 
izontal displacement eigenfunctions, and p the local density. 
It follows that g-modes will have higher kinetic energies than 
p-modes, as they propagate in denser regions of the star. 
This is clearly seen by comparing Fig.[2]and the correspond- 
ing plot for p-modes, Fig. 1101 

Following the asymptot i c the o ry of nonrad i al os- 
cillations (see e.g. | Tassoull Il980l ; ISmevers et al.l 1 19951 : 
ISmevers fc Movall2007l ). modes of the same degree I and 
consecutive radial order k would have a uniform period sep- 
aration, AP, given by: 



AP = Pi 



fc+i,( 



Pk,i = 2tt 



\N\ 



dr 



(2) 



where N is the BV frequency. Thus, mode trapping effects 
are revealed for certain modes as deviations with respect to 
the mean period spacing. This can be seen from the data 
points plotted in red in Fig. [3] 

The radial and horizontal displacement eigenfunctions 
(yi = £ r /r and yp = ( a 2 r/g) * £h/r, respectively as denned 
by iDziembowskil Il97ll ) of representative modes with mini- 
mum (g88), normal (g85) and maximum (g83) total kinetic 
energy are shown in Fig. [4] The vertical line at log q ~ —2.23 
marks the location of the maximum composition gradient 
in the He/H transition region. It is evident that the 388 
mode has lower amplitude below the He/H interface, so it 
is trapped in the envelope. As the region below the He/H 
transition is very dense, it has a high weight in the total ki- 
netic energy, and hence trapped modes have minimum val- 
ues of the kinetic energy. We recall that, for high-radial or- 
der ^-modes in this model, significant energy interchange is 
produced in the region below the He/H transition, which is 
mostly damping (fig. 8, Paper I), resulting in trapped modes 
being less damped and having maximum growth rate values. 
On the other hand, the g83 mode has highest amplitude in 
the damping region, and it is referred to as a confined mode, 
which is less likely to be excited, as is reflected in its mini- 
mum value of the growth rate. 

Previous adiabatic studies for sdBs (|Charpinet et al.1 

|2000| ) postulated qualitatively, based on the growth rate de- 
pendence of the kinetic energy, r\ a W/ Ekin (where W is the 
work integral denned as the total energy balance over one 
period of oscillation), that when a mode had lower ampli- 
tude in the damping regions, the mode would be less damped 
and would be more likely to be excited. Indeed, we confirm 
this hypothesis for high-radial order sdO gravity modes in 
our non-adiabatic study. 

We also draw the attention to the fact that the radial 
eigenfunction for the trapped mode has a node above the 
He/H transition, while the confined mode has a node below 
it, and the normal mode has a node just about on the in- 
terface. For the horizontal eigenfunction, although there are 
nodes at similar distances above and below the He/H tran- 
sition, for the trapped (confined) mode, the closest node to 
the transition is below (above) it, in good agreement which 
what is found for sdBs and white dwarfs. However, the larger 
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Figure 4. Radial (loft) and horizontal (right) displacement eigcnfunctions for trapped (388), normal (g85) and confined ($83) g- modes 
of the original model. The vertical line marks the location of the He/H transition. The insets show a blow-up of this region. 
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Figure 5. Weight function for 388, gS5 and gS3, a trapped, normal and confined mode, respectively for the original model (left) and 
modes with equivalent frequencies for the model without He/H chemical transition (right). Inset showing a blow-up of the He/H transition 
region. 



distances of the nodes to the trapping region in our model 
translate into less efficiency of the trapping/confinement 
mechanism. 

In our model, the smaller differences in the kinetic en- 
ergy between trapped and confined modes comparing to sdB 
or white dwarf modes, are due to the fact that the mode 
trapping interface (the He/H transition, see below) is lo- 
cated much deeper in the envelope (log q ~ —2) than for 
sdBs (logq ~ — 4. in | Charpinet et al.|[2000l ) or white dwarfs 
(logg ~ —10, in lBrassard et alj|l992h . As a result, both the 
trapped and confined region, weighted by the density, p (see 
Eq. [TJ , contribute largely to the kinetic energy. 

Also from ICharpinet et al.l (|2000l ) and iBrassard et al.l 

(| 19921 ) , trapped modes in sdBs and white dwarfs have sharp 
kinetic energy minima, while maxima, corresponding to con- 
fined modes, are somewhat wider, comprising up to three 
modes and indicating that confinement processes are not 
so effective as mode trapping. In our model, both the min- 
ima and the maxima in the kinetic energy have a certain 
width, revealing that both the mode trapping and confine- 
ment processes are not so sharp as in sdBs and white dwarfs. 



The explanation lies in the wider chemical transition regions 
in the sdO model (see Fig. [T] right) in co mparison to those 
for sd Bs and whit e dwarfs (see e.g. fig. 3 of lCharpinet et all 
l200d and fig. 1 of lBrassard et alj|l992h . 

A clearer view of the local contribution of each re- 
gion of the star to the formation of a particular mode 
is given through the so-called weight function, defined in 
ICharpinet et all (|2000l ) as: 



F(^ r ,P',$';r) = 



£n 2 + {F " 



TiPp 



+ 



P' 
TiP 



TV 2 

+ 6- — 



pr 2 (3) 



where P' and $' are the Eulerian perturbations of the 
pressure and potential, respectively and Ti is an adiabatic 
exponent. 

A plot of the weight function for the same trapped (g88), 
normal (#85) and confined mode (#83) (Fig. [S] left), each of 
them normalised individually, shows that the wide C-O/Hc 
chemical transition region has an important contribution for 
all of them, as expected for high-radial order g-modes, which 
mainly propagate in the deep interior of the star. The main 
difference in the weight function for the different modes is 
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Figure 6. Left: Logarithm of the kinetic energy vs. radial order for low-radial order g-modes of the original (upper circles and dotted 
line) and perturbed models without the He/H (upper solid line) and C-O/Hc (upper crosses and dashed line) transition regions. Lower 
symbols plot the corresponding growth rates whose values can be read on the right axis. Right: Period spacing for low-radial order 
g-modes with the same degree i and consecutive radial orders for the original and perturbed He/H and C-O/He models (with same 
symbols as in the left pannel). 




Figure 7. Original (dotted) Brunt- Vaisala frequency (left) and sound speed (right) and modifications cancelling out the He/H (solid) 
and sharp peak of the C-O/He (dashed) chemical transition regions. 



the higher contribution of the He/H transition region to the 
formation of the trapped mode, which also has a higher con- 
tribution of the outer zones. This becomes even clearer when 
we eliminate in the model the He/H chemical transition re- 
gion (Fig. [5] right, explanation below). 

A blow up of Fig. [2] for the low-radial order g^modes (k 
up to 30) is shown in the left panel of Fig. [6] (circles and dot- 
ted lines) , and the period differences in this zone in the right 
panel. We still notice the existence of certain modes alter- 
nating maximum and minimum values of the kinetic energy, 
until about k — 10. The period differences between modes 
also show deviations with respect to the mean value, a signa- 
ture of mode trapping. However, contrary to high-radial or- 
der g-modes, this trapping/confinement phenomenon, does 
not have any influence on the ability to excite the modes, as 
they remain highly stable. This is explained as, for g-modes 
of low-radial order, the energy interchanged below the He/H 
transition is negligible compared to the energy interchanged 
at the Z-bump region, located at —10.0 < logg < —9.5. 



There, the energy contribution for both trapped and con- 
fined modes is always damping (see fig. 8 in Paper I). So, 
even if the mechanical effect of the mode trapping is main- 
tained, at the frequencies in which low-radial order g-modes 
occur, the significant energy interchange is produced in dif- 
ferent regions of the star, therefore the trapping has no in- 
fluence in the growth rate. 

3.1 Perturbation analysis 

We have investigated the effects of cancelling out the He/H 
and the sharp peak of the C-O/He chemical transition 
regions of the Brunt- Vaisala frequency and sound speed 
(Fig. [7} in this sdO model on mode trapping and the ten- 
dency of modes to instability. To do this, we described the 

3 These new perturbed models will be referred to as pcrtHcH and 
pertCO, and their symbol codes in the plots are solid lines, and 
crosses and dashed lines respectively. 
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Figure 8. Differential work for g88 mode of the original (dotted) 
model and g87 mode of the pertHeH (solid) model. The inset 
zooms in on the He/H chemical transition. See text for details. 
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Figure 9. Growth rates of the original model (circles), the same 
model with the differential work integrated only over log 9 ^ 
—2.24 (crosses) and the pertHeH model (solid line). 



perturbation of the mass, pressure, sound speed squared (c 2 ) 
and BV squared (-/V 2 ) as a function of the perturbation of 
the density, gravity and adiabatic exponent. From these re- 
lations we derived an expression for the perturbation of the 
density as a function of the perturbations introduced for 
BV and the sound speed (SN 2 , 8c 2 ). Finally, we wrote the 
perturbation of the dynamical variables of the model as a 
function of the already calculated perturbations (a complete 
description of the mathematical formulation is given in Ap- 
pendix A). These new perturbed models underwent again 
the non-adiabatic analysis. 



3.1.1 Cancelling out the He/H transition 

When we cancel out the He/H chemical transition region, we 
get a smooth kinetic energy for all the modes (Fig. [2] upper 
solid line). Trapped modes vanish, revealing that the He/H 
transition is the main factor responsible for their occurrence. 
The plot of the weight function (Fig. [5] right) reveals that 



Figure 10. Logarithm of the total kinetic energy vs. radial order 
for p-modes of the original (upper circle and dotted line) and 
perturbed models without the He/H (upper solid line) and C- 
O/Hc (upper crosses and dashed line) transition regions. Lower 
symbols plot the corresponding growth rates whose values can 
be read on the right axis. Filled circles mark p4 and p9 modes 
showing maximum growth rate values. See Section 5 for more 
details. 



eto650t45 1 = 2 p — modes 




k (1 = 2) 

Figure 11. Large frequency spacing vs. radial order for p-modes 
with the same degree i and consecutive radial orders for the orig- 
inal, and perturbed models. 



the He/H region, once cancelled, no longer has any signif- 
icant influence on mode formation. We also note that the 
growth rate no longer oscillates with frequency (Fig. [21 lower 
solid line), as expected from the uniform kinetic energy; and 
we also get a uniform period separation (Fig. [3j solid line). 
However, perturbed modes equivalent to trapped modes in 
the original model achieved higher values of the growth rate, 
while their kinetic energy was higher, which seems to be at 
odds with the growth rate dependence on the kinetic energy. 

The local minimum in the damping energy for model 
8 at logg = —2.23 (Fig. [8] dotted line), which is absent 
in model pertHeH (Fig. [S] solid line), is responsible for the 
overall lower values of the growth rate in the original model. 
Thus, when we calculate the work function for the original 
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model integrating between log q ^ —2.24 and the surface, we 
obtain similar growth rate values as those of the pertHeH 
model (Fig. EJ). 

In the case of low radial order (/-modes (Fig. H3 solid 
lines) some subtle mode trapping seems still to be present. 
This may be due to remnant mode trapping effects caused by 
the C-O/He transition, in the same way as described below 
for low radial order p-modes. 

3.1.2 Cancelling out the C-O/He transition 

When we cancel out the sharp peak of the C-O/He chem- 
ical transition, keeping the He/H peak, we still obtain the 
original pattern of trapped modes, with total kinetic en- 
ergy, growth rate and period differences similar to those of 
the original model (Fig. [2] Fig. [3] crosses and dashed lines) . 
Thus, this chemical transition plays no role in high radial or- 
der g-mode trapping. However, for lower order modes (Fig. [6] 
left panel, crosses and dashed lines), it may have some in- 
fluence, which is revealed in certain alterations to the total 
kinetic energy of each mode. 



and confined mode (Fig. 1131 right) shows that significant 
energy interchange is only produced at the Z-bump; and 
there, driving is higher for the trapped mode. As it is the 
case for low-radial order (/-modes, a maximum in kinetic en- 
ergy does not translate into a minimum growth rate value for 
the mode. This is explained by the maximum kinetic energy 
being due to a higher amplitude of the eigenfunctions in the 
region logq > —0.8, which for p-modes does not have signif- 
icant influence on driving. Neither do the higher amplitudes 
of the eigenfunctions of the trapped mode at the driving 
region lead to a maximum growth rate value, as consecu- 
tive higher radial order modes have even higher amplitudes. 
We conclude that the influence of the work function prevails 
over the influence of the kinetic energy in the computation 
of the growth rate. 

Pressure modes with higher radial orders, not shown 
in Fig. 1101 display kinetic energy monotonically increasing, 
an effect of the nodes accumulating in the surface limit im- 
posed by the bo u ndary conditions, as it is well described in 
ICharpinet et al.l (|2000h . The trapping effect caused by the 
C-O/He transition is no longer produced, as is expected for 
high-radial order pure p-modes propagating only in the ex- 
ternal layers of the star. 



4 PRESSURE MODES 

The logarithm of the total kinetic energy, and corresponding 
growth rates, for the p-mode spectrum are shown in Fie. 1101 
The depicted radial orders correspond to frequencies rang- 
ing from about 6 to 25 mHz. There are certain modes which 
show values higher than the mean kinetic energy in the orig- 
inal and perturbed models. We investigate if this could also 
be an effect of mode trapping. 

Following agai n the a sympt o tic theory of non- 
radial oscillations ([Tassoull Il980l ; ISmevers et al.l 1 19951 ; 
ISmevers fc Moval l2007f T p-modes of the same degree I and 
consecutive radial order k, will have constant frequency sep- 
aration, Av, given by : 



P 

Fip 



1/2 



dr 



Thus, the mode trapping signature for p-modes would 
show up in the deviations in the large frequency spacing for 
the original and the perturbed He/H model (Fig. 1110 . 

We have plotted the radial and horizontal displacement 
eigenfunctions for the original model p-modes with maxi- 
mum (p9), intermediate (p8) and minimum (p7) total kinetic 
energy (Fig. I12|l . We found that all modes show higher am- 
plitude in the envelope, as expected for p-modes. We found 
that the p9 mode shows the highest amplitude in the outer 
envelope for the 3 modes under comparison. However, ul- 
timately what accounts for the higher total kinetic energy 
is the higher amplitude of the eigenfunctions in the inner- 
most layers of the star, where the mode is trapped due to 
the pinching effect of the C-O/He interface on the eigen- 
functions. This can be seen in Fig. [13] (left) that shows the 
local contribution to the kinetic energy of the trapped, nor- 
mal and confined modes, p9, p8, p7, respectively, normalised 
to the maximum of p7. Although not shown for the sake of 
clarity, the maximum kinetic energy value for the p9 trapped 
mode reaches one order of magnitude higher. 

The plot of the differential work for a trapped, normal 



4.1 Perturbation analysis 

When we cancel the He/H, or the sharp peak of the C- 
O/He chemical transitions in the Brunt- Vaisala frequency, 
and also the corresponding bumps in the sound speed at 
these locations, we still find a kinetic energy pattern, and an 
oscillating profile of the growth rate, both similar to those 
of the original model (Fig. 1101 solid and dashed lines, re- 
spectively). Therefore, neither the He/H, nor the C-O/He 
chemical transition, nor the sound speed at these locations, 
seem to have any significant influence on the kinetic energy 
or on the tendency to driving of the p-modes. 

This result was the expected for the growth rate, for 
reasons given above, but not for the kinetic energy, as we 
presumed the modification of the C-O/He transition would 
show in the kinetic energy profile, due to changes in the am- 
plitude of the eigenfunctions. However, eigenfunction pro- 
files similar to those of Fig. [T2] are retained for trapped, 
normal and confined modes of the perturbed models. The 
reason may be due to that, in fact, the C-O/He transition 
still exists (see Fig. [7] left) as we only cancelled out its steep 
peak, giving a softer profile. 



5 SUMMARY 

We conclude for the g-mode spectrum that: 

• High-radial order (/-modes are trapped in the envelope 
due to the He/H chemical transition. This provides a weak 
selection mechanism, in the way that trapped modes oscil- 
late with lower amplitudes in the innermost damping region. 
Thus, trapped modes show lower total kinetic energy and 
higher values of the growth rate. 

• Low to intermediate radial order g-modes (mixed 
modes) are also trapped in the envelope by the He/H transi- 
tion. However, this has no significant effects on the driving, 
as the innermost damping region, at the frequencies at which 
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Figure 12. Radial (left) and horizontal (right) displacement eigenfunctions for a trapped, normal and confined p-mode of the original 
modcl:p9, p8 and p7, respectively. Dashed-dotted vertical lines delimit the convcctivc region where maximum opacity takes place. Dashed 
vertical lines in the insets mark the maxima of the C-O/He (the deepest) and Hc/H transition regions. 




Figure 13. Left: Local kinetic energy for trapped, normal and confined modes: p9, p8, p7, respectively. Right: Differential work and 
opacity (long dashed line) for the same modes. The maximum of the Z-bump is located at logg ~ —9.5. 



low g-modes occur, has no significant weight in the energy 
interchanged by the mode with its surroundings. 

• Cancelling out the He/H transition region cancels the 
mode trapping effects of high order gravity modes, and the 
growth rate loses its oscillatory behaviour with frequency, as 
expected. However, modes with local maxima in the kinetic 
energy and similar interchanged energy at the driving re- 
gion than their non-perturbed counterparts, achieve higher 
growth rate values. This effect is due to the presence of a 
small but strong damping region in the original model at the 
location of He/H chemical transition, which has an overall 
damping effect in the high-radial order g-modes of the origi- 
nal model. For low radial order gravity modes some remnant 
trapping still exists which may be attributed to the C-0 /He 
transition. 

• When we cancel out the sharp peak of the C-O/He tran- 
sition region, keeping the He/H peak, mode trapping is still 
present for high order g-modes, and also the oscillating be- 
haviour of the growth rate. Thus, the C-O/He transition 
region has no significant role in their trapping. However, it 
may play a subtle role in low order g-mode trapping (seen 



in certain alterations to the total kinetic energy) , due to the 
persistence of the transition. Trapping though, has no effect 
in their driving. 

We also have to bear in mind that although the Brunt- 
Vaisala, frequency is expected to have a higher weight in the 
mode trapping, the contribution of the sound speed has not 
been disentangled in the analysis and may be responsible in 
part for certain aspects of trapping. 

We conclude for the p-mode spectrum that: 

• Low to intermediate radial order p-modes (mixed 
modes) are found to be trapped by the C/O-He transition 
in the innermost region. However, this has no significant ef- 
fects on the driving as the maximum energy interchanged is 
produced in the envelope of the star. 

• Cancelling out the He/H or the sharp peak of the C- 
O/He transition has no significant effect neither in the mode 
trapping, nor in the growth rate profile, which remain es- 
sentially unaltered. This is interpreted as trapping effects 
produced by the remaining piece of the C-O/He transition. 
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6 DISCUSSION AND CONCLUSIONS 

A sdO model presented an oscillating profile of the growth 
rate with frequency, which lead us to analyse it in more de- 
tail: it was found that gravity modes suffered mode trapping 
effects caused by the He/H chemical transition region which 
made high-order ^-trapped modes more likely to be driven. 
Mode trapping was also found for low-radial order p-modes, 
caused by the deep C-O/He chemical interface. However, 
no obvious correlation could be made with its possible non- 
adiabatic effects. 

We sought to investigate the non-adiabatic effects of 
mode trapping through the influence of the kinetic energy 
on the growth rate beh aviour. Up to now, the existing adi- 
abatic studies for sdBs (|Charpinet et al]|2000l ) qualitatively 
proposed that, as a trapped mode had lower amplitude in 
the damping regions, the mode would be less damped and 
would be more likely to be excited. Indeed, this is what we 
found in our non-adiabatic study for trapped in the envelope 
high-radial order p-modes: they oscillate with lower ampli- 
tudes below the He/H transition, a damping region, hence 
they have higher values of the growth rate than its confined 
counterparts. In this way, mode trapping acts a weak selec- 
tion mechanism rendering trapped modes more likely to be 
driven. 

Meanwhile, low- to intermediate-radial order p-modes, 
which show a similar behaviour of the growth rate, were 
found to be trapped by the C-O/He transition. However, as 
the significant energy interchange at p-mode frequencies is 
produced much higher in the envelope, mode trapping under 
the C-O/He interface does not play a significant role in the 
potential to destabilise modes. 

Therefore, we conclude that mode trapping is a poten- 
tial aid for destabilizing modes when (1) modes are trapped 
in such a way that they have lower amplitude in a damp- 
ing region and (2) this region plays a significant role in the 
energy interchanged by the mode with its surroundings. 

Finally, a tentative exercise of a theoretical mode iden- 
tification of the radial order of the observed modes of 
JI600+0748 was done in the terms of asymptotic analysis, 
based on the similar mean frequency spacing of our model 
8 and the observed frequencies. We associated the latter 
with those showing highest growth rate in model 8, and 
gaps in the observed frequencies to modes with minimum 
growth rate. This resulted in the observed modes associ- 
ated with high radial-order p-modes, in contradiction with 
low-order low-degree mode identification by Fontaine et al. 
(2008). However, we consider this possible application of 
mode trapping to explain the observed sdO pulsation spec- 
tra, a path that is worth exploring. For this exercise to be 
fruitful, proper theoretical models reproducing the observed 
physical parameters of J1600+0748 should be used. 



Spanish Ministry of Education and Science. RO is supported 
by the Research Council of Leuven University through grant 
GOA/2003/04. The research of JM is supported in part by 
a grant from the Mount Cuba Astronomical Foundation. 



APPENDIX A 

We describe here the modifications imposed on the dynam- 
ical variables of the model through perturbations of Brunt- 
Vaisala, frequency (5N 2 ) and the sound speed (5c 2 ). The 
strategy that we follow is first to describe the perturbation 
of physical variables of mass, pressure, N 2 and c 2 as a func- 
tion of the perturbation of the density (Sp), gravity (Sg) and 
adiabatic exponent (<5Fi). Then, we derive an expression for 
the perturbation of the density as a function of 5N 2 , 5c 2 . 
Finally, we will describe the perturbation of the dynamical 
variables of the model as a function of the already calculated 
perturbations. 

We begin perturbing the equations of stellar structure: 
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We perturb the sound speed: 
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Using the equations we have derived, we operate each 
term in the equation above as: 
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Replacing these three terms in Eq. [TT]we obtain: 
Sg 



g 5c gdpSp g dSp 
c 2 c 2 p dr p p dr 



Replacing JV 2 by Eq. \W\ we obtain the integro- 
differential equation in Sp: 
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Substituting Eq. [10] in the last term, replacing Sg by 
Eq. [S]and using Eq. Owe find: 



P ( 2g 2 9 dp\ 8g 
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p dr 
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Numerical tests indicate that this last term is negligi- 
ble. Thus, we are left with a linear first order differential 
equation of the type y + a(x)y + b(x) = with solution: 
y = JcMKW^ with c{x) = g/ »(«)*._ ThuSj solving for 

Sp we get: 
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In this way we have the perturbed density Sp as a func- 
tion of the perturbation of Brunt- Vaisala, <5iV 2 and the sound 
speed Sc 2 . 

Finally, we perturb five independent dynamical vari- 
ables of the twenty four input variables which describe our 
equilibrium models. These input variables are the same 
whic h describe the equilibrium m odels in adipls adiabatic 
code (|Christensen-Da lseaard 2008) and are convenient w hen 
the equations are formulated as in iDziembowskil (|l97ll ). as 
it is our case: 
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From Eq. |9] we have the perturbation of the adiabatic 
exponent: 
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